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In this lecture, we will focus on configurations where the angular frequency of the light is close to 
some transition frequency of the medium. In particular, we will start with a brief outline of how 
the non-resonant susceptibilities may be modified in such a way that weakly resonant interactions 
can be taken into account. Having formulated the susceptibilities at weakly resonant interaction, 
we will proceed with formulating a non-perturbative approach of calculation of the polarization 
density of the medium. For the two-level system, this results in the Bloch equations governing 
resonant interaction between light and matter. 

The outline for this lecture is: 

• Singularities of the non-resonant susceptibilities 

• Alternatives to perturbation analysis of the polarization density 

• Relaxation of the medium 

• The two-level system and the Bloch equation 

• The resulting polarization density of the medium at resonance 

Singularities of non-resonant susceptibilities 

In the theory described so far in this course, all interactions have for simplicity been considered 
as non-resonant. The explicit forms of the susceptibilities, in terms of the electric dipole moments 
and transition frequencies of the molecules, have been obtained in lecture six, of the forms 

X^i (-w; uj) ~ ^ba_ + { similar [B. & C. (4.58)] 

i Iba - U) 

X { "U-^;co u w 2 ) ~ , r ^7 -v + {similar terms}, [B.&C. (4.63)] 

x^ 7 (-^;^,c 2 ,c 3 ) ~ m <^£c£k + {simi i ar tcrms}; 

[B.&C. (4.64)] 



To recapitulate, these forms have all been derived under the assumption that the Hamiltonian 
(which is the general operator which describes the state of the system) consist only of a thermal 
equilibrium part and an interaction part (in the electric dipolar approximation), of the form 

H = H + Hi(t)- 

This is a form which clearly does not contain any term related to relaxation effects of the medium, 
that is to say, it does not contain any term describing any energy flow into thermal heat. As long 
as we consider the interaction part of the Hamiltonian to be sufficiently strong compared to any 
relaxation effect of the medium, this is a valid approximation. 

However, the problem with the non-resonant forms of the susceptibilities clearly comes into 
light when we consider an angular frequency of the light that is close to a transition frequency of 
the system, since for the first order susceptibility, 

Xm2(-^; w ) -> °o, whenu)~>tt ba , 
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or for the second order susceptibility, 

X^api-^'^' w z) ->■ 00, when uj 1 + uj 2 -*■ ^ba or u 2 -> rj ca - 

This clearly non-physical behaviour is a consequence of that the denominators of the rational 
expressions for the susceptibilities have singularities at the resonances, and the aim with this 
lecture is to show how these singularities can be removed. 

Modification of the Hamiltonian for resonant interaction 

Whenever we have to consider relaxation effects of the medium, as in the case of resonant interac- 
tions, the Hamiltonian should be modified to 

H = H + H I (t)+H R , (1) 

where, as previously, Hq is the Hamiltonian in the absence of external forces, H\(t) = —Q a E a (r, t) 
is the interaction Hamiltonian (here taken in the Schrodinger picture, as described in lecture four), 
being linear in the applied electric field of the light, and where the new term £/r describes the 
various relaxation processes that brings the system into the thermal equilibrium whenever ex- 
ternal forces are absent. The state of the system (atom, molecule, or general ensemble) is then 
conveniently described by the density operator formalism, from which we can obtain macroscopi- 
cally observable parameters of the medium, such as the electric polarization density (as frequently 
encountered in this course), the magnetization of the medium, current densities, etc. 

The form (1) of the Hamiltonian is now to be analysed by means of the equation of motion of 
the density operator p, 

iH^=Hp-pH=[H,p}, (2) 

and depending on the setup, this equation may be solved by means of perturbation analysis (for 
non-resonant and weakly resonant interactions), or by means of non-perturbative approaches, such 
as the Bloch equations (for strongly resonant interactions). 

Phenomenological representation of relaxation processes 

In many cases, the relaxation process of the medium towards thermal equilibrium can be described 

by 

[Hn,p] = -ihf(p-po), 

where po is the thermal equilibrium density operator of the system. The here phenomenologically 
introduced operator T describes the relaxation of the medium, and can can be considered as being 
independent of the interaction Hamiltonian. Here the operator T has the physical dimension of an 
angular frequency, and its matrix elements can be considered as giving the time constants of decay 
for various states of the system. 

Perturbation analysis of weakly resonant interactions 

Before entering the formalism of the Bloch equations for strongly resonant interactions, we will 
outline the weakly resonant interactions in a perturbative analysis for the susceptibilities, as pre- 
viously developed in lectures three, four, and five. 

By taking the perturbation series for the density operator as 

p(t) = Jh^ + faQ) + p 2 {t) +...+ p n (t) +..., 

~[E(t)]° -[^(t)] 1 ~[E(t)] 2 ~[E(t)]» 
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as we previously did for the strictly non-resonant case, one obtains the system of equations 



dt 
1 dt 



[H ,p 2 (t)} + [H I (t),p 1 (t)]-iHfp 2 (t), 



ih ^df L = $0,Pn(t)] + [Hl(t),Pn-l(t)]-ihtp n (t), 



As in the non-resonant case, one may here start with solving for the zeroth order term p , with all 
other terms obtained by consecutively solving the equations of order j = 1, 2, . . . , n, in that order. 

Proceeding in exactly the same path as for the non-resonant case, solving for the density 
operator in the interaction picture and expressing the various terms of the electric polarization 
density in terms of the corresponding traces 

oo oc 



P,(r,t) = ][>("> (r,i) = I]TTr[p„(t)CU 



V 

n=0 n=0 

one obtains the linear, first order susceptibility of the form 



,.(1) I , ,. , ,\ _ iVe _ ( n \ ( r ab r ba , r ab r ba \ 



Similarly, the second order susceptibility for weakly resonant interaction is obtained as 

(2) I \ 



e a h 2 2! ^ be + -iT ac )(n ab + uj a -iT ab ) 

ab be ca ab be' ca 



(Q ac + LU2- «r ac )(O bc + LO a - iT bc ) (tt ba +LU 2 - iT ba )(^bc + W<r - i^bc 



r r a n 
1 ab 1 be 1 ca 



(tt ba +LU 2 - iT ha ){^ca + ^ - iT ca ) . 

In these expressions for the susceptibilities, the singularities at resonance are removed, and the 
spectral properties of the absolute values of the susceptibilities are described by regular Lorenzian 
line shapes. 

The values of the matrix elements T mn are in many cases difficult to derive from a theoretical 
basis; however, they are often straightforward to obtain by regular curve-fitting and regression 
analysis of experimental data. 

As seen from the expressions for the susceptibilities above, we still have a boosting of them close 
to resonance (resonant enhancement). However, the values of the susceptibilities reach a plateau 
at exact resonance, with maximum values determined by the magnitudes of the involved matrix 
elements T mn of the relaxation operator. 



3 



Nonlinear Optics 5A5513 (2003) 



Lecture notes 11 



Validity of perturbation analysis of the polarization density 

Strictly speaking, the perturbative approach is only to be considered as for an infinite series 
expansion. For a limited number of terms, the perturbative approach is only an approximative 
method, which though for many cases is sufficient. 

The perturbation series, in the form that we have encountered it in this course, defines a power 
series in the applied electric field of the light, and as long as the lower order terms are dominant in 
the expansion, we may safely neglect the higher order ones. Whenever we encounter strong fields, 
however, we may run into trouble with the series expansion, in particular if we are in a resonant 
optical regime, with a boosting effect of the polarization density of the medium. (This boosting 
effect can be seen as the equivalent to the close-to-resonance behaviour of the mechanical spring 
model under influence of externally driving forces.) 

As an illustration to this source of failure of the model in the presence of strong electrical fields, 
we may consider another, more simple example of series expansions, namely the Taylor expansion 
of the function sin(x) around x « 0, as shown in Fig. 1. 

2 1 1 1 1 1 1 1 



1.5 




sin(x) « Pz{x) = x 



Figure 1. Approximations to f(x) = sin(x) by means of power series expansions of various degrees. 

In analogy to the susceptibility formalism, we may consider x to have the role of the electric 
field (the variable which we make the power expansion in terms of), and sin(x) to have the role of 
the polarization density or the density operator (simply the function we wish to analyze) . For low 
numerical values of x, up to about x « 1, the sin(x) function is well described by keeping only the 
first two terms of the expansion, corresponding to a power expansion up to and including order 
three, 

X 3 

3! ■ 

For higher values of x, say up to about x « 2, the expansion is still following the exact function to 
a good approximation if we include also the third term, corresponding to a power expansion up to 
and including order five, 

x 3 x 5 

m + 5T 

This necessity of including higher and higher order terms goes on as we increase the value of x, 
and we can from the graph also see that the breakdown at a certain level of approximation causes 
severe difference between the approximate and exact curves. In particular, if one wish to calculate 
the value of the function sin(x) for small x, it might be a good idea to apply the series expansion. 



sin(x) « pz(x) = x 
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For greater values of x, say x w 10, the series expansion approach is, however, a bad idea, and an 
efficient evaluation of sin(a;) requires another approach. 

As a matter of fact, the same arguments hold for the more complex case of the series expansion 
of the density operator 1 , for which we for high intensities (high electrical field strengths) must 
include higher order terms as well. 

However, we have seen that even in the non-resonant case, we may encounter great algebraic 
complexity even in low order nonlinear terms, and since the problem of formulating a proper 
polarization density is expanding more or less exponentially with the order of the nonlinearity, the 
usefulness of the susceptibility formalism eventually breaks down. The solution to this problem is 
to identify the relevant transitions of the ensemble, and to solve the equation of motion (2) exactly 
instead (or at least within other levels of approximation which do not rely on the perturbative 
foundation of the susceptibility formalism). 

The two-level system 

In many cases, the interaction between light and matter can be reduced to that of a two-level 
system, consisting of only two energy eigenstates \a) and \b). The equation of motion of the 
density operator is generally given by Eq. (2) as 



with 

H = H + H(t) + H R . 

For the two-level system, the equation of motion can be expressed in terms of the matrix elements 
of the density operator as 

ih^ = [H ,p] aa + [Hl(t),p]aa + [Hn,p]aa, (3a) 
if%d W = [ ^f>U + [Hi(t),p} ab + [H R ,p] ab , (3b) 
ih d E^_ = [H ,p] bb + [Hi(t),p] bb + [H R ,p\ bb , (3c) 



where the fourth equation for p ba was omitted, since the solution for this element immediately 
follows from 

Pba = Pah- 
Terms involving the thermal equilibrium Hamiltonian 

The system of Eqs. (3) is the starting point for derivation of the so-called Bloch equations. Starting 
with the thermal-equilibrium part of the commutators in the right-hand sides of Eqs. (3), we have 
for the diagonal elements 

[H ,p\aa = (a\H p\a) - (a\pH \a) 

= ^{a\Ho\k){k\p\a) - ^{a\p\j) (j|£o|a> 

= E a /7 aa — PaaM-a 
= 

= [Hoip]bb, 



1 We may recall that the series expansion of the density operator is the very origin of the expansion 
of the polarization density of the medium in terms of the electric field, and hence also the very 
foundation for the whole susceptibility formalism as described in this course. 
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[H ,p} ab = (a\H p\b)-(a\pH \b} 

= ]T H£4fc)(fc|p|6) - 5>|p|j>£!£o|6>, 

fc =E 5 afe j =V. } Sjb 

= EaPab - Pab^b 

= -(E fc -E a ) Pafc 

= -Ml ba p ab 

Terms involving the interaction Hamiltonian 

For the commutators in the right-hand sides of Eqs. (3) involving the interaction Hamiltonian, we 
similarly have for the diagonal elements 

[Hi(t),p] aa = (a\(-er a E a (r,t))p\a) - (a\p(-ef a E a (r,t))\a) 

= -e J B a (r,i)|^(a|f Q |fc}<fc|p|a}-^(a| / 3|i)0-|f a |a} 
^ k j 

= -eE a (r,t)^r2 a Paa + r" b p ba - p aa xZ a - PabCj 

= - e ( r abPba - r% a p ab )E a (r, t) 
= -[H I (t)J] bb , 

and for the off-diagonal elements 

[#i(*),pU - (a\(-ef a E a {r,t))p\b) - (a\p(-ef a E a (r,t))\b) 

= - eJ E Q (r,i){E( a l^l A: )( fc l/ 3 l 5 )-E( a l^)0'l^l 6 ) 

^ k j 

= -eE a (r, i) jOab + r" b p bb - p aa r* b - p o6 rg, j 

= -er" b E a (r,t)(p bb - p aa ) - e(r* a - r% b )E a {v,t)p ab 

= {Optical Stark shift : SE k = -er% k E a (r, t), k = a, 6} 

= -er% b E a {r,t)(p bb - p aa ) + (SE a - SE b )p ab . 

Terms involving relaxation processes 

For the commutators describing relaxation processes, the diagonal elements are given as 

[H R ,p\ aa = -ih(p aa - p a {a))/T a , 
[H R ,p] bb = -ih(p bb - p {b))/T b , 

where T a and T b are the decay rates towards the thermal equilibrium at respective level, and where 
Po(a) and po(b) are the thermal equilibrium values of p aa and p bb , respectively (i. e. the thermal 
equilibrium population densities of the respective level). The off-diagonal elements are similarly 
given as 

[H R ,p] ab = -ihp ab /T 2 , 
[H K ,p] ba = -ihp ba /T 2 . 

A common approximation is to consider the two states \a) and \b) to be sufficiently similar in 
order to approximate their lifetimes as equal, i. e. T a « T b « Ti, where T\ for historical reasons 
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is denoted as the longitudinal relaxation time. For the same historical reason, the relaxation time 
T 2 is denoted as the transverse relaxation time. 2 

As the above matrix elements of the commutators involving the various terms of the Hamilto- 
nian are inserted into the right-hand sides of Eqs. (3), one obtains the following system of equations 
for the matrix elements of the density operator, 

ih ~ir = ~ e ( r " bPba ~ r baPab)E a (r, t) - ih(p aa - po(a))/T a , (4a) 
ih ~^F = ~ hriba P ab ~ er ab E a{r,t)(pbb - Paa) + (SE a - SE b )p ab - ihp ab /T 2 , (4b) 
ih< ^r = e ( r * bPha ~ r baPab)E a {v, t) - ih(p bb - Po(b))/T b . (4c) 

(The system of equations (4) corresponds to Butcher and Cotter's Eqs. (6.35).) So far, the applied 
electric field of the light is allowed to be of arbitrary form. However, in order to simplify the 
following analysis, we will assume the light to be linearly polarized and quasimonochromatic, of 
the form 

E a {r,t) = Re[E%{t)eM-iwt)]. 

We will in addition assume the slowly varying temporal envelope E"(t) to be real- valued, and we 
will also neglect the optical Stark shifts <5E a and 5E b . In the absence of strong static magnetic 
fields, we may also assume the matrix elements er" b to be real-valued. When these assumptions 
and approximations are applied to the equations of motion (4), one obtains 

i{,Pba - Pab)P(t) cos(wt) - (p aa - p (a))/T a , (5a) 



dt 
dpab 
dt 

dpbb 
dt 



i^baPab + lf3(t) COs(wt)(p bb - Paa) ~ Pab/T 2 , (5b) 

-i{pba - Pab)P(t) cos(wt) - (p bb - po(b))/T b , (5c) 



where the Rabi frequency /3(t), defined in terms of the spatial envelope of the electrical field and 
the transition dipole moment as 



/?(*) = er« b EZ(t)/h = er ab ■ E w (i)//i, 



was introduced. 



The rotating-wave approximation 

In the middle equation of the system (5), we have a time-derivative of p a b in the left-hand side, 
while we in the right-hand side have a term itlbaPab- Seen as the homogeneous part of a linear 
differential equation, this suggests that we may further simplify the equations of motion by taking 
a new variable p^ b according to the variable substitution 

Pab = Pab CXp[i(n ba - A)t], (6) 

where A = ilb a - w is the detuning of the angular frequency of the light from the transition 
frequency tt ba = (E& - E a )/H. 



2 For a deeper discusssion and explanation of the various mechanisms involved in relaxation, 
see for example Charles P. Slichter, Principles of Magnetic Resonance (Springer- Verlag, Berlin, 
1978), available at KTHB. This reference is not mentioned in Butcher and Cotters book, but it is 
a very good text on relaxation phenomena and how to incorporate them into a density-functional 
description of interaction between light and matter. 
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By inserting Eq. (6) into Eqs. (5), keeping in mind that p ba = P ab , one obtains the system 
dp a 



dt 
dt 

dpbb 
dt 



i{p\ l a exp[-i(Q 6a - A)t] - p l ab exp[i(n ba - A)t])/?(i) cos(wt) - ( Paa - p (a))/T a , (6a) 
iAp% b + 0(t) cos(wi) cxp[-i(n ba - A)t]( Pbb -p aa )- p^ b /T 2 , (6b) 
-i(pg, cxp[-i{n ba - A)t] - p% b cxp[i(n ba - A)t])[3(t) cos(wi) - ( Pbb - P o(b))/T b , (6c) 



The idea with the rotating-wave approximation is now to separate out rapidly oscillating terms of 
angular frequencies lo + Q ba and — (u) + f2& a ), and neglect these terms, compared with more slowly 
varying terms. The motivation for this approximation is that whenever high-frequency components 
appear in the equations of motions, the high-frequency terms will when integrated contain large 
denominators, and will hence be minor in comparison with terms with a slow variation. In some 
sense we can also see this as a temporal averaging procedure, where rapidly oscillating terms 
average to zero rapidly compared to slowly varying (or constant) components. 

For example, in Eq. (6b), the product of the cos(wi) and the exponential function is approxi- 
mated as 

cos(wi) exp[— i(fl ba — A)t] = ^-[cxp(iu;t) + cxp(— iuit)] cxp[—i (fl ba — A) t] 

2 s „ / 

— OJ 

= ^[l + exp(-i2wt)] -> \, 
while in Eqs. (6a) and (6c), the same argument gives 

exp[i(fl ba — A)t) cos(cji) = -[cxp(iuit) + cxp(-iujt)} cxp[— i (Q ba — A) t] 

2 s „ / 

— LJ 

= ^[exp(i2wi) + l]-i 

By applying this rotating-wave approximation, the equations of motion (6) hence take the form 

^ = \(fit P*b)m (Paa Po(a))/T a , (7a) 

^ = iA P % b + l -f3(t)(Pbb Paa) P° b /T 2 , (7b) 
^ = ~\(Pba - PaoW) (Pbb P0(b))/T b . (7C) 

In this final form, before entering the Bloch vector description of the interaction, these equations 
correspond to Butcher and Cotter's Eqs. (6.41). 

The Bloch equations 

Assuming the two states \a) and |6) to be sufficiently similar in order to approximate T a « T b « T l7 
where T x is the longitudinal relaxation time, and by taking new variables (u,v,w) according to 

o i ft 

U = P ba + Pab, 

■ i n n \ 

V = nPba-Pab)> 
W = Pbb- Paa, 

the equations of motion (7) are cast in the Bloch equations 

— = -Av-u/T 2 , (8a) 
-=Au + f3(t)w-v/T 2 , (8b) 
d ^ = -P(t)v-(w-w )/T 1 . (8c) 
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In these equations, the introduced variable w describes the population inversion of the two-level 
system, while u and v are related to the dispersive and absorptive components of the polarization 
density of the medium. In the Bloch equations above, w = pa (b) — po (a) is the thermal equilibrium 
inversion of the system with no optical field applied. 

The resulting electric polarization density of the medium 

The so far developed theory of the density matrix under resonant interaction can now be applied 
to the calculation of the electric polarization density of the medium, consisting of N identical 
molecules per unit volume, as 

P ll (T,t) = N(ef lt ) 

= NTr\per li ] 

= NJ2 ik\per»\k) 

k—a.b 

= N H E mmer»\k) 

k—a,b j—a,b 

= N E {{k\p\a)(a\ef^\k) + (k\p\b){b\er,\k)} 

k—a.b 

= N{(a\p\a)(a\ef^\a) + (b\p\a)(a\ef^\b) + (a\p\b) <6|ef» + {b\p\b) {b\ef ^b)} 
= N(p ba er^ b + p ab er^ a ) 

— {Make use of p ab = (u + iv) cxp(iu>t) — pl a } 

= N[(u — iv) cxp(—iujt)er^ b + (u + iv) exp(iu>t)er ba \. 

The temporal envelope P£ of the polarization density, throughout this course as well as in Butcher 
and Cotter's book, is taken as 

P"(r,t)=Re[P^exp(-iwt)], 

and by identifying this expression with the right-hand side of the result above, we hence finally 
have obtained the polarization density in terms of the Bloch parameters (u, v,w) as 

P£{T,t) = NerZ b (u-iv). 

This expression for the temporal envelope of the polarization density is exactly in the same mode 
of description as the one as previously used in the susceptibility theory, as in the wave equations 
developed in lecture eight. The only difference is that now we instead consider the polarization 
density as given by a non-perturbative analysis. Taken together with the Maxwell's equations (or 
the proper wave equation for the envelopes of the fields), the Bloch equations are known as the 
Maxwell-Bloch equations. 
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